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ABSTRACT 

The supranova model of gamma-ray bursts (GRBs) , in which the GRB event is preceded by a supernova 
(SN) explosion by a few months to years, has recently gained support from Fe line detections in the X- 
ray afterglows of some GRBs. A crucial ingredient of this model yet to be studied is the fast-rotating 
pulsar that should be active during the time interval between the SN and the GRB, driving a powerful 
wind and a luminous plerionic nebula. We discuss some observational consequences of this precursor 
plerion, which should provide important tests for the supranova model: 1) the fragmentation of the 
outlying SN ejecta material by the plerion and its implications for Fe line emission; and 2) the effect 
of inverse Compton cooling and emission in the GRB external shock due to the plerion radiation field. 
The plerion-induced inverse Compton emission can dominate in the GeV-TeV energy range during the 
afterglow, being detectable by GLAST from redshifts z <^ 1.5 and by AGILE from redshifts z <J 0.5, and 
distinguishable from self-Compton emission by its spectrum and light curve. The prospects for direct 
detection and identification of the precursor plerion emission are also briefiy considered. 

Subject headings: gamma rays: bursts, pulsars: general, stars: neutron, supernova remnants, radiation 
mechanisms: non-thermal, line: formation 



1. INTRODUCTION 

In the currently standard interpretation of gamma-ray 
bursts (GRBs), the central engine gives rise to a highly 
relativistic outflow, the 'fireball'. Although the pres- 
ence of relativistic outflows in GRBs has been amply 
demonstrated by multiwavelength observations of after- 
glows (Piran 1999, van Paradijs, Kouveliotou & Wijers 
2000, Meszaros 2001), the nature of the central engine 
itself is still a great mystery. One model of the central 
engine that has recently gained attention is the 'supra- 
nova' model of Vietri & Stella (1998, hereafter VS98), 
in which a massive star ends up in a supernova (SN), 
but the subsequent black hole formation and GRB event 
is delayed by some time tp, typically expected to be a 
few months to years. The SN is assumed to leave be- 
hind a rotationally-supported, 'supramassive' neutron star 
(SMNS; hereafter also referred to simply as 'pulsar') which 
then slowly shrinks by shedding angular momentum via a 
magnetospheric wind (and/or gravitational radiation; see 
§5). After a spin-down time tp when roughly half its ini- 
tial angular momentum has been lost, the configuration 
becomes unstable and collapses to a black hole, possibly 
with a surrounding disk, leading to the GRB proper.^ The 
model's major advantage lies in the potential realization 
of a very baryon-clean pre-GRB environment (mandatory 
for generating sufficiently relativistic fireballs): first the 
SN ejects the majority of the outlying mass of the pro- 
genitor star, and second the SMNS driven-wind effectively 



sweeps up remaining baryonic matter in the vicinity of the 
central object. Support for this model comes from recent 
detections of strong Fe emission features in the X-ray af- 
terglow spectra of some GRBs (Piro et al. 1999, Antonelli 
et al. 2001, Yoshida et al. 2001), particularly that of 
GRB991216 (Piro et al. 2000), as weU as a transient Fe 
absorption feature in the prompt emission of GRB990705 
(Amati et al. 2000), indicating surprisingly large amounts 
of Fe-rich material existing nearby, yet relatively removed 
from the GRB site. Observational constraints on its lo- 
cation, quantity, density and velocity are compatible with 
the pre-ejected shell SN remnant (SNR) in the supranova 
model (Lazzati et al. 1999, Vietri et al. 2001, Lazzati et al. 
2001), but may be difficult to accomodate in other mod- 
els (see however Rees & Meszaros 2000, Meszaros & Rees 
2001, Bottcher & Fryer 2001, McLaughhn et al. 2002). 

The precursor pulsar wind should be extremely pow- 
erful. Before collapsing, the SMNS must inevitably ex- 
pel a substantial fraction of its rotational energy Ep as a 
magnetically-driven wind. 



E„ 



GjM^uj 
2^ 



2.4 X lO^^erg jo.eM^tOi, 



(1) 



where M = SMsMq and uj = W^ujiS'^ are the SMNS's 
mass and angular velocity, j — 0.6j'o.6 is its angular mo- 
mentum in units of GM^ / c, and the numbers are typical 
SMNS model values (Salgado et al. 1994, Cook, Shapiro 
& Teukolsky 1994, VS98). The spin-down time (i.e. the 

^The numerical simulations of SMNS collapse to a black hole studied by Shibata, Baumgarte & Shapiro (2000) do not support the formation 
of a debris disk that may power the GRB by accretion or external electromagnetic torques. However, the GRB may also be triggered through 
the winding-up of strong magnetic fields due to differential rotation of the SMNS during the course of its collapse (e.g. Kluzniak & Ruderman 
1998). 
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SN-GRB delay time) and corresponding wind luminosity 
can be estimated from the magnetic dipole formula (Pacini 
1967, VS98) 

*^ = B^R^LO^ - ^^'^^^^ jo.6MlR-^,,w^^B-l^ , (2) 

and 

= Ep/t, = 7.1 X lO^^erg Rl^^uisl,^) 

where i?, = 15i?*,i5km is a typical SMNS equatorial ra- 
dius. The surface magnetic field = 10'^'^i?*,i3G is un- 
constrained from model calculations and can be consid- 
ered a free parameter. Equivalently, we may take tp as 
the free parameter, and vary Lp = Ep/tp accordingly with 
Ep fixed as in eq.l; the fiducial value we choose below is 
tp ^ 120days (implying Lp ^ 2.3 x lO^^erg s~^) so as to 
be consistent with observations of GRB991216 (see §3.3). 
During tp, Lp is expected to be relatively constant, and the 
wind should energize a plerionic nebula in the pre-GRB 
surroundings, a more compact yet much more luminous 
version of the Crab nebula. 

The consequences of such a precursor plerion in the 
supranova scenario has not been considered previously, 
and this paper addresses some important dynamical and 
radiative effects it may induce, each providing important 
observational diagnostics for the supranova model. We dis- 
cuss the acceleration and fragmentation of the SN ejecta 
material by the plerion-SNR interaction and its implica- 
tions for Fe line emission in §2, and inverse Compton 
scattering of the ambient plerion radiation field in the 
GRB external shock and the resulting high-energy after- 
glow emission in §3. A brief consideration of the direct 
detection and identification of the precursor plerion emis- 
sion is given in §4. We will assume a fiat lambda cosmology 
with = 0.3, flA = 0.7 and Hq = TOkms^^Mpc'^ 

After this paper was submitted, we became aware 
of the work of Konigl and Granot (2002, hereafter 
KG02), who also investigated the properties of GRB af- 
terglows occurring inside plerionic nebulae in the context 
of the supranova model. They stressed the advantages 
of such a picture in realizing the relatively high mag- 
netic fields and electron injection efficiencies required in 
the GRB blastwave to explain observed afterglows, and 
demonstrated this by constructing magnetohydrodynami- 
cal (MHD) models for the plerion under various assump- 
tions. Here we are interested in some characteristic ob- 
servable effects caused by the precursor plerionic activity 
which are peculiar to the supranova model and through 
which the model may be tested by futm-e observations. 
The emphasis of this paper will be on the relevant radiative 
processes and their observational implications; we choose 
to keep the discussion of the plerion dynamics relatively 
simple and leave more detailed modeling of this aspect to 
future studies. A study following our approach but with 
a more realistic treatment of the plerion has recently been 
carried out by Guetta & Granot (2002). 

2. DYNAMICAL EFFECTS OF THE PLERION ON THE 
SUPERNOVA REMNANT 

An exemplary case of the dynamical interaction between 
a plerion and a SNR can be seen in the well-studied Crab 



nebula, which is known to be accelerating and fragment- 
ing the surrounding SNR, resulting in the prominent op- 
tical filaments instead of a clear shell (e.g. Davidson & 
Fesen 1985, Hester et al. 1996). In this respect, the pow- 
erful supranova plerion can be even more effective than 
the Crab. We model the plerion in a simple way fol- 
lowing Pacini & Salvati (1973, hereafter PS73; see also 
Bandiera, Pacini & Salvati 1984, Chevalier & Reynolds 
1984, Chevalier 2000), considering a homogeneous, spheri- 
cal bubble into which energy is injected at a constant rate 
Lp = W^^LpAQerg s"\ a fraction = 0.5^b,o.5 going 
into magnetic field and the rest = 1 — = 0.5^e,o.5 
into relativistic electrons. The electrons here are mostly 
radiative (see §3), and we neglect their pressure for sim- 
plicity. As the plerionic bubble initially plows through 
the expanding core and envelope of the progenitor star, it 
should accelerate the swept-up ejecta material (Chevalier 
1977, Chevalier & Fransson 1992, KG02). If a fair fraction 
of the total plerion energy Ep ~ lO'^'^erg can be conveyed 
to SN ejecta of several M©, their final attained velocities 
could reach Vs ~ 0.05 — 0.1c, consistent with those inferred 
from the observed width of the Fe line in GRB991216 (Piro 
et al. 2000). 

As in the Crab, the SN ejecta should not entirely re- 
main as a spherical shell during this acceleration phase 
due to Rayleigh- Taylor (RT) instabilities operating at the 
plerion-SNR interface. The growth timescale of the RT in- 
stability on a spatial scale R is tj^x ^ {R/RY^'^^ and since 
R ~ AttR'^pb/Ms, pb = Bp/Sn being the plerion magnetic 
pressure, its ratio to the expansion timescale texp ^ R/vg 
is 

tflT (2M^y^' f AE, y/' 

te.p^\BlR^j ^\3^bEpJ 

~ 0.39Elg(CB,o.5ip,46)-'/'iro (4) 

(Bandiera, Pacini & Salvati 1983), where Ms, Vg and 
Eg = Msvl/2 = 10^^i?s,5ierg are the SNR ejecta mass, 
velocity and kinetic energy, and tio is time after the SN 
in units of 10 days. In the case of the Crab, tpT/texp > 1 
and the RT instability is currently setting in, consistent 
with other lines of evidence (Hester et al. 1996, Sankrit 
& Hester 1997). However, it can potentially develop much 
more rapidly for the supranova plerion. 

The final outcome of the RT instability is not easy to 
predict, as it will depend on the instability's non-linear 
behavior as well as other processes such as radiative cool- 
ing and thermal conduction. We refer to the work of Jun 
(1998) as a guideline, who carried out 2-dimensional nu- 
merical simulations of the plerion-SNR interaction and as- 
sociated RT instabilities, choosing Crab-like parameters 
for Lp and t and neglecting cooling. His results demon- 
strate that the RT instability leads to dramatic effects at 
late times, strongly disrupting the SNR shell and trans- 
forming it into pronounced RT "fingers" which carry the 
majority of the SN ejecta mass and kinetic energy. To be 
identified with the Crab's optical filaments, these fingers 
protrude inward from the outer ejecta shell and end in 
dense "heads" with slower expansion velocities than the 
plerion, which can eventually become decoupled from the 
shell. In Jun's simulations, the density of the heads was 
on average ^ 10 times that of the rest of the SN ejecta. 
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and the radial extent of the fingers was ^ 0.2 times that 
of the sheU radius. However, as mentioned in his diseus- 
sion, the inclusion of radiative cooling should increase the 
density contrast of the heads, and both this and the ra- 
dial extent of the fingers should grow with time. These 
quantities could then be much larger for the supranova 
plerion, since the relative timescale of the instability here 
might be much shorter (eq.4). Although a firm conclusion 
warrants high resolution simulations of the relevant condi- 
tions, it is not implausible that the plerion can effectively 
shred a substantinal portion of the SN ejecta into con- 
densed fragments and engulf them during its lifetime tp. 
Note that Arons (2002), in a somewhat different context, 
has recently discussed the issue of SNR fragmentation by 
a powerful pulsar wind in greater depth. 

If this is indeed the case (see §5 for comments when 
this is not), there are important implications for the Fe 
line emission in afterglows. The simplest explanation for 
the lines is Fe Ka miiltiple recombination radiation from 
Fe-rich matter photoionized by the X-ray afterglow con- 
tinuum on timescales of days (Lazzati et al. 1999, Weth 
et al. 2000, Paerels et al. 2000). Consistency with obser- 
vations require that the material to be illuminated is very 
dense and possesses a large covering factor, and yet that 
the GRB line of sight is devoid of any such matter so that 
the GRB blastwave will not be decelerated too quickly 
(Lazzati et al. 1999). The filamentary fragmentation of 
the SNR may naturally account for such a geometry. 

The condensation of the SNR matter may also en- 
hance the line emissivity. The emission rate of line pho- 
tons through recombination when a sufficient photoion- 
izing flux irradiates material with electron density Ug = 
IG^ne.gcm"^ and temperature T = lO^T^K containing an 
Fe mass Mpe = O.lMFe.o.iM© is 

Nline = A^FeArec ^ 1-5 X lO^'^S'^ MFe,0.ine,9Ti^^^ (5) 

where A^pc is the total number of Fe atoms and tree — 
127s Tj^'^n'^g is the recombination time for H-like Fe 
(Paerels et al. 2000, Vietri et al. 2001). By the time 
of the GRB at tp = 120tp,i2odays (see §3), a SNR of 
Ms = IOMs^iqMq expanding at Vs = O.lvgfi.ic will have 
reached a radius Rg — 3.1 x 10^^t;s,o.iip,i2ocni, and if dis- 
tributed in a spherical shell of width AR^ ^ O.lRg, its 
electron density would be 

ne,s ^ 5.3 X 10^MsAa{vs,o.itp,i2oy^cm~^. (6) 

Typically observed line fluxes can be obtained for 
Mpcfi.i ;> 4 (Eq.5), but only if the entire spherical shell 
is irradiated and emits line photons efficiently. The GRB 
blast wave must then be isotropic and no beaming is al- 
lowed (Piro et al. 2000; see however KG02), which is in- 
congruous with recent evidence to the contrary (e.g. Frail 
et al. 2001, Panaitescu & Kumar 2002). However, if 
the SNR cjccta can be condensed into filaments of nnich 
higher density, tree is shortened, and irradiating a smaller 
fraction of material may suffice to produce the observed 
lines. For comparison, the density of the the Crab's opti- 
cal filaments averaged over the volume of the whole neb- 
ula is n 10cm~^, whereas the actual density of the 

^EC processes have been extensively discussed in relation to blc 
shocks, see e.g. Bottcher & Dermer 1999. 



individual filaments is n ^ lO'^cm^"^ (Davidson & Fesen 
1985). Although detailed numerical modeling is neces- 
sary for quantitative predictions, we may speculate that 
the density contrast could be even higher for the supra- 
nova plerion due to the instability being more developed 
(sec above), and beaming fractions of order ^ 0.01 may 
be entirely compatible with the observed Fe lines. Note 
that in explaining the transient Fe absorption feature in 
GRB990705, Lazzati et al. (2001) have inferred the pres- 
ence of clumpy SNR ejecta material with density con- 
trasts of order 100-1000 surrounding the burst site (see 
also Bottcher, Dermer & Fryer 2002). The observed line 
energy in GRB991216 can also be consistent with beam- 
ing if identified with blueshifted emission from He-like Fe 
(for which the rest frame energy is 6.7 keV, and tree is 
similar to H-like Fe at somewhat lower T) in material ap- 
proaching toward the observer with Vg ^ 0.1c (Ballantyne 
& Ramirez-Ruiz 2001). 

We mention that even when circumstances conducive 
to Fe line emission cannot be realized, the ac;tion of the 
plerion for tp ^ lO'^s may expel much of the baryons sur- 
rounding the GRB site out to radii ^ lO^^cm and clear the 
path for the ensuing fireball, which is a great advantage of 
the supranova model. 

3. EXTERNAL COMPTON COOLING AND EMISSION DUE 
TO THE PLERION RADIATION FIELD 

Another important implication of the precursor plerion 
in the supranova model is that the GRB should be ex- 
ploding into a 'radiation-rich' environment, i.e. into the 
luminous radiation field of the plerion, which can consti- 
tute a large fraction of the pulsar wind luminosity. Such 
'external' photons impinging into the relativistically ex- 
panding GRB shock from outside would be seen highly 
Doppler-boosted in the shock comoving frame, and act as 
efficient seeds for inverse Compton scattering (referred to 
as external Comptonization, or EC), with important ob- 
servational consequences for the afterglow emission. ^ 

Since our focus here is on the possible manisfestations 
of the EC emission process in the supranova model, we 
do not attempt a detailed physical modeling of the the 
plerion dynamics, employing instead a simplified descrip- 
tion for illustrative purposes. As discussed by KG02, such 
modeling would entail various uncertainties and requires a 
dedicated study by itself, which is out of our current scope. 
Once a suitable dynamical model of the plerion can be 
constructed, the formulation of this paper may be used to 
make more reliable predictions of the afterglow emission. 
Alternatively, if an afterglow EC emission component is 
actually observed in the future, we may interpret the data 
with our formulation so that the unknown properties of 
the supranova plerion can be probed. 

3.1. Plerion Emission 

We continue the discussion of the plerion as a roughly 
homogeneous, spherical bubble with constant energy in- 
jection. Following the arguments of the preceding section, 
it is postulated that the plerion can first penetrate effec- 
tively through the expanding SN ejecta, shredding and 
entraining the ejecta material within a relatively short 

;; see e.g. Sikora 1997. For previous work on EC processes in GRB 
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time. Thereafter, the plerion can be considered to ex- 
pand into the; ambient medium, here taken to be uniform 
for brevity. The consequences of the dense, entrained 
baryonic clumps for the global dynamics of the plerion 
is unclear, but lacking a realistic model, we do not delve 
into this complication and make simplified assumptions 
regarding its evolution. If the SNR filaments are expand- 
ing at velocities of order Vs ^ 0.1c (as implied by the 
observed Fe line width in GRB991216), the expansion 
velocity of the plerionic bubble Vp may be several times 
larger than this, judging from the observed expansion ve- 
locities in the Crab as well as the numerical simulations 
of Jun (1998). As long as the ejecta material is not dy- 
namically well-coupled to the plerion as a whole, we may 
also estimate Vp from the relativistic blastwave solution 
of Blandford & McKee (1976) in the case of steady, adia- 
batic energy injection into a uniform medium, which gives 



surface of the plerion is 

JvP 



(7) 



Synchrotron self-absorption will cause an effective mini- 
mum cutoff in the spectrum at a frequency 



ttx/Btt 



2nmpC 



2.4 X 10 



3^/(s + l)eifei?p 
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where f{s -I- 1) is a factor of order unity weakly depending 
on s (e.g. Rybicki & Lightman 1979), and 



K, 



{kLp/p,c-f/H-''^ 



1/4,-1/2 



p,120i 



where 



o-TcCBipln(7M/7m) 



(9) 



is the normalization factor of the cooled electron distri- 
bution (PS73). The expressions in the last lines of eqs.8 
and 9 are for our fiducial choice of s = 2. As the flux 
falls steeply below u^, the overall plerion spectrum is that 
of eq.7 between and i^^, with constant luminosity per 
logarithmic frequency interval, and the total emitted lu- 
minosity is of order £,eLp. '"^ Note that some; quantities 
regarding the plerion contain subscripts or superscripts p 
in order to avoid confusion with those for the GRB blast- 
wave. 

We mention that the plerion should also drive a forward 
shock into the ambient medium, whereby the associated 
synchrotron radiation may give rise to an additional emis- 
sion component. Since the magnetic field and the electron 
injection efficiency in the shocked external medium are 
expected to be lower than that inside the plerion, such 
emission should be of subdominant luminosity appearing 
at frequencies below u^. EC scattering of this compo- 
nent by GRB blastwave electrons is unimportant for the 
observed afterglow, as it is likely to be masked by the 
synchrotron and SSC radiation occurring in the pertinent 
frequency range (see below); however the emission at radio 
frequencies may play a role in the direct detectability of 
the plerion before the GRB (§4). A proper account of this 
plerion forward shock emission mandates a more detailed 
description of the plerion expansion, so in this study we 
choose to neglect this component and restrict ourselves to 
emission from the plerion interior. 

3.2. GRB Afterglow Emission 

At time tp after the SN, the pulsar collapses and the 
GRB goes off, sending a fireball and relativistic blastwave 
into the plerion. We may distinguish two different situa- 
tions for how the deceleration of the GRB blastwave, and 
hence the afterglow, initiates. As in conventional discus- 
sions of afterglows, the deceleration may begin when it has 
swept up enough outlying baryonic material, which may 
occur if the (probably coUimated) fireball in our line of 
sight happens to strike the denser parts of the entrained 

^This formula strictly applies only in the ultrarelativistic limit (Fp 1, Vp ^ c); in the nonrelativistic limit (Tp ~ 1, <S c), the 

deceleration goes as Vp oc tp and Rp oc tp^'' . 

*With our fiducial choice of €s = €e = 0.5, synchrotron-self- Compton cooling is at most comparable to synchrotron cooling, so this is 
neglected for simplicity. 

^Any flux above which may arise from self-Compton emission is irrelevant for electrons in the GRB shock due to Klein-Nishina effects; 
see below. 



Tp — (1 — Vp/c?) is the expansion bulk Lorentz factor, 

ni = pi/nipC^ = ni^ocm"*^ is the ambient density and k is 
a numerical factor of order unity. ^ A mildly relativistic 
expansion is then inferred for ages tp ^ 10^ — lO^days and 
the values of Lp and ni considered here. We assume, at 
least for the conditions of our interest, that the plerion has 
been expanding at a nearly constant velocity of Vp ~ c, so 
that its radius Rp ~ dp ~ 3.1 x 10-^''tp_i2ocm. 

The luminosity and spectrum emitted by the plerion 
may be evaluated following PS73. The magnetic field 
inside the plerion (taken to be isotropically tangled) is 
^, o cr<ci/2 +^i^/2^ with the value 

of Ep fixed as in eq.l. Relativistic electrons/positrons 
(hereafter simply 'c;k;ctrons') are injected into the plerion 
with a power of ^eEp, having an isotropic, power-law en- 
ergy distribution dn/djP oc (7^)"'* in the Lorentz factor 
range < < 7^,^ . We fiducially take ^5 = = 0.5, 
s = 2 and 7^^ = 1 , parameters which lead to consistent fits 
when similar models are applied to the multiwavelength 
spectra of the Crab and other known plerionic nebulae 
(PS73, Bandiera et al. 1984, Amato et al. 2000, Chevalier 
2000). The maximum Lorentz factor 7^^ is set by equat- 
ing the synchrotron cooling time t^y ~ Gnmec/arBp'- 
with the acceleration time, here assumed to be t^^^ 
2Trjpmec/eBp (consistent with Fermi acceleration at rel- 
ativistic shocks), so that 7^ ~ {Se/arBpy^"^ and the 
corresponding maximum synchrotron emission frequency 
z/^ ~ Se'^/nmeCaT ^ 1.2 x lO^^Hz. ^ The distribution 
should develop a break at the Lorentz factor 7P where t^ 



Rp/vp 



t 



pi 



equals the adiabatic expansion time t^^ 
leading to 7P ^ GTrniec/arBptp. Our fiducial parameters 
give 7^ ^ 6, implying that injected electrons of all but 
the lowest energies arc radiative and cool within t^^; the 
cooled distribution has index s + 1 and emission spectrum 
of energy index Up ^ s/2 = 1. The specific flux at the 



5 



SN cjecta fragments inside tlie plerion. Alternatively, the 
deceleration may result from the blastwave acciimulating 
enough inertia in relativistic particles and magnetic field 
of the pulsar wind material, as discussed by KG02. To de- 
scribe the afterglow in this case, KG02 showed that stan- 
dard formulae (e.g. Sari, Piran & Narayan 1998) can still 
be used if one replaces parameters such as the ambient 
medium density and electron and magnetic energy den- 
sity fractions by equivalent quantities related to the ple- 
rion properties (see below). In considering either case, we 
will assume for simplicity that the decelerating medium 
is roughly uniform, at least within the timescales we con- 
sider. (Note that this is not a bad approximation for most 
of the MHD plerion models discussed by KG02.) 

The standard expressions for the radius r and bulk 
Lorentz factor F of the shocked material in an adiabatic, 
spherical blastwave decelerating self-similarly in a uni- 
form medium are r(t) = [\2Ect/ Aimm,pC?{l + z)]^/^ ^ 
3.6 X 10i^cm(^52/n)i/'*[td/(l + z)]^/'^ and T{t) = [3^(1 + 
zf /2f>&T^nmpcH^)Y/^ ~ 5.9(i;52/n)V8[t^/(l + 2)]-3/8, re- 
spectively (e.g. Meszaros & Rees 1997, Piran 1999). Here 
E = 10'^^i?52erg is the blastwave energy, ncm^'^ the exter- 
nal medium density, t = t^day the observer time elapsed 
after the GRB, z is the GRB redshift, and we have adopted 
the kinematic relation t = r{l + z)/4r'^c. The decel- 
eration starts at radius r^ec = (3i5/47rnmpC^rQ)^/3 ~ 

at observer time tdec — 

rdec(l + z)/4Tlc ~ 2.4s(l + z){E52/n)^/^T-%l where 
To = 300ro,3oo is the initial bulk Lorentz factor. 

To describe the time-dependent, multiwavelength after- 
glow spectrum from the blastwave, we follow standard 
discussions of the synchrotron emission (e.g. Meszaros 
& Rees 1997, Sari, Piran & Narayan 1998, Wijers & 
Galama 1999), and extend it to include cooling and emis- 
sion by the EC process in addition to the synchrotron-self- 
Compton (SSC) process (Panaitescu & Kumar 2000; Sari 
& Esin 2001, hereafter SEOl; Zhang & Meszaros 2001, 
hereafter ZMOl). It is assumed that constant fractions 
cb = 0.01eB,-2 and Se = 0.1ee,-i of the postshock en- 
ergy are imparted to magnetic field and relativistic elec- 
trons, respectively. The comoving magnetic field is B(t) = 

{32Tre BUmpy/^Tc ~ 0.23G€]i^_^Eli\^/»[td/il + z)]-^/». 
Electrons are accelerated in the shock to a power-law 
distribution dn/d'y oc in the Lorentz factor range 

7m < 7 < 7m- We assume p > 2, our standard choice 
being p = 2.5. The minimum Lorentz factor 7^ is given 

by 



rripp- 



ruep-l 

360ee,-i[3(p-2)/(p-l)] 

x{E,2/n)'/'[ta/{l+z)]-^/' 



(10) 



(note that [3{p - 2)/(p - 1)] = 1 when p = 2.5). The 
electrons radiatively cool by the combination of the syn- 
chrotron, SSC and EC processes, the timescales of which 



are t' 



Yt',.y and tL 



Xt' 



sy 



respectively, the total cooling time being t'^ — [{l/t'gy) 
(1/Cc) + (l/C)]-' = C(l + Y + ^)"'- The maximum 



Lorentz factor 7m is determined by balancing t'^ with the 
acceleration time t^^^ 2iTjmec/eB, 



lM{t) 



3e 



1/2 



arBl + Y + X ^ 
~ 0.97 X 10^(1 + Y + X)-^'^e-]L\ 

x£;-//%-3/16[t,/(l + ^)]3/16 (11) 

(ZMOl). The cooling Lorentz factor 7c is where t'^ equals 
the comoving adiabatic expansion time ~ r/cT ~ tV, 



1 



l + Y + X gtTBH 



2.9 X IQ^{1 + Y + X)-^eg]_^ 
xE^,'/\-^/^[Ul{l + z)Y'\ 



(12) 



The observed characteristic synchrotron emission fre- 
quency for an electron of Lorentz factor 7 is = 
(4/3)r(3/47r)(eB/mec)72(l -|- z)'^ (Wijers & Galama 
1999), so those for each of 7^, 7m and 70 are respectively 

Vrn{t) ~ 1.0 X 10i2Hze]/'_2e^,_i[3(p- 2)/(p- 1)]2 



UM{t) ~ 0.72 X 1023Hz(1 + Y + X)-'^ 

xEli\-''h-/'\l + z)-'/\ (14) 



and 



Vc{t) - 6.4 X 10^''^Hz(l + r + X)-'^eg^!3^ 
^El^"n-H-^l\x^zr-l\ 



(15) 



At early times {t < tt), jm > 7c, and all electrons 
cool within t'^^ (fast cooling regime). The synchrotron 
flux then peaks at i^,- and mainly consists of three 
power-law segments: f„ = (;^/fc'c)^/^/i/,max for ly < Uc, 
= {v/vc)'^^"^ fv,ma^ for Vc < I' < I'm and = 

(i^m/i^c)"^^^(i^/i^m)~^^Vi^,max for < V < I'M- ^ At latC 

times (t > it), 7m < 7c, and only electrons of 7 > 7c cool 
within (slow cooling regime). Then the spectrum peaks 
at Vm, and again has three parts: = {v /vmY^^ fv,max for 

V < Vm: fu = (i^/i^m)"'''"^'^^/i.,max for i^,„ < < l^c and 

U = [vclvmY'^'^^^^^'^iy IvcY'^^'^ for < V < VM. 
The peak synchrotron flux /^,max for either fast or slow 
cooling is 

/.,ma. = CrBr3(l + z)/i^i 

~ 2.9mJye^/i2£;52ni/2z)-2g(l + ^), (16) 

where C is a numerical constant and = 10^8-DL,28cm is 
the luminosity distance to the GRB. The transition from 
fast to slow cooling (7c = 7^) occurs at the time 



tt 



13s(l + Y + Xfe%^_,el_MP ' 2)/(p - 1)]' 

xE52n{l + z).{17) 



^Synchrotron self-absorption should cause an additional spectral brecik (or possibly two breaks for fast cooling; Granot, Piran & Sari 2000) 
at the lowest frequencies, but is ignored here as it is observationally irrelevant for the high energy afterglow emission, our main concern. 
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To account for the effects of EC and SSC cooling, we 
have introduced X and Y, which are in general time- 
dependent. The EC parameter X equals the ratio of co- 
moving frame energy density in plerion radiation field u'p 
to that in magnetic field u'^. If the blastwave decelerates 
inside the plerion (r < Rp), which is the case mainly con- 
sidered below, we may presume that the radiation field is 
nearly uniform and isotropic, obtaining 



X 



p2 £,eLp 



21n(7L/7^) 



where v'^j^ ~ rrieC? /h'y is the frequency above which 
Klein-Nishina (KN) effects suppress Compton scattering, 
depending on 7. Observe here that the factors of in 
u'p and u'^ cancel out so that aside from the weak t- 
dependence through v'^^, X is constant in time. Gen- 
erally v'^ > but we do not make large errors in 
taking min{f^, i/^j^} u'^ in the logarithmic factor 
for simplicity. For a blastwave decelerating outside the 
plerion (r > -Rp), the radiation field is anisotropic and 
falls with r, and X differs from eq.l8 by a factor of 
~ {Rp/rf{l -2n + V)/3 where ^JL = ^1 - {Rp/rf 
(which applies for F » 1; see e.g. Sikora et al. 1996, 
Inoue & Takahara 1997), decreasing with time. 

The SSC parameter Y can be evaluated as in SEOl, 
y ~ <yh'B ~ m'JWBO^+y+X)] ~ rieJ[eB{l+Y+X% 
where u'^ and u'^y are the comoving frame energy densi- 
ties in relativistic electrons and synchrotron radiation, and 
77 ^ min{l, (7c/7m)^^^} is the fractional energy radiated 
by electrons. In the fast cooling phase, 77 ~ 1 and we get 
the expression 



y-\U^^^Hi^xY-i-x]. 



(19) 



In the slow cooling phase, Y decreases with time but only 
very slowly for typical values of p (SEOl), so we simplify 
by assuming eq.l9 at all times. 

The SSC spectrum has a similar shape to the syn- 
chrotron spectrum, and we approximate it as comprising 
broken power-laws (SEOl, ZMOl), with characteristic fre- 
quencies at 



sc 2 



1.2 X I0''mel_,e'i%[2,{p - 2)/{p - l)f 
xE%\-^'H-/'\l + zf/^ (20) 



and 



~ 5.4 X 10^'^Hz(l + Y + Xy^e^'^^ 

•KEl^l^n-''l%^l\\ + z)-^l\ (21) 
The maximum SSC frequency should occur at 
vkn ^ F7MmeC^(l z)~^ 

~ 0.71 X 1028Hz(1 + y + X)-i/2e-^/2 

Xi;5^/^«n-5/16t-3/16(i + ^)-13/16^ (22) 



where KN effects completely cut off the spectrum. (Gen- 
erally, Vkn ^ vff ~ Jm^m-) For fast coohng, the SSC 
flux Z^*^ peaks at I'f", the spectrum being oc v^^^ for 
p < v^^ , oc for vf-^ < V < v'^ and oc i/"^/^ for 

< v < Vkn', for slow coohng, the peak is at v^ , 
and the spectrum is oc v^^^ for u < v^ , oc z/~(p~i)/2 



for 



,sc 



< V < v^'^ and oc v~'p/'^ for v^'^ < v < 



Vkn- 



The peak SSC flux /,f tiqual to the peak synchrotron 
/i/,max multiplied by the optical depth Tg of the shocked 
material, 



/AeBT^nmpC^, (18) 



sc 

max 



arnr 



2.3 X 10-'m3ye]i'EU'n'/'DlU'J\l + 



The EC spectrum reflects the shapes of both the elec- 
tron distribution and the plerion spectrum. Considering 
deceleration inside the plerion, and again adopting the bro- 
ken power-law approximation, spectral breaks arise at the 
frequencies where electrons of 7^ and 7c upscatter plerion 
photons of pP, 



,EC 



jlT'vP{i + zy 

~ 1.1 X lO^^Hze? 



nl/2 - 



1/2. 



^,_i[3(p-2)/(p-l)p 
'/'(l + ^)i/2 (24) 



and 



0.71 X 10^^Hz(l + Y + Xye} 



xE-,'/\-^/X 



1/2 



(1 + ^) 



-1/2 



(25) 



where the fiducial numerical value for v^ in eq.8 has been 
used. Note the factor F^ which results from two Lorentz 
transformations, one into the comoving frame from the 
observer frame and another vice-versa. The fast cool- 
ing EC spectrum has a peak flux at v^'-', is oc for 

V < v^^ from the behavior of the Compton scattering 
cross section, oc v^'^/'^ for vf-^ < v < v^ and oc for 
v^ < V < Vkn, mirroring the flat plerion spectrum. The 
KN limit discussed above also applies here for the max;- 
imum EC frequency, v^^ being irrelevant. Likewise, the 
slow cooling EC spectrum peaks at v^ , and is oc for 

V < vE^, oc for v^^ < V < v^^ and oc v''^ 
for vf^ < V < VKN- The peak EC flux /^^^x niay be 
obtained analogously to the SSC case, but using the re- 
lation /^f j^ax ~ '''ef'v'p max ^hc comoving frame, where 
f'v'p max ~ ^fvi- max the the peak plerion flux emitted at 
v'^ = v'P ^ Fi'^. Accounting for Lorentz transformations 
to and from the comoving frame, 



f^c- ^ p2 p _!_(-! I ,^ 

J !^,max ' i/.max j-j2 ' 

3 87ri?2ln(7^/7P,)i5| 



(1 + ^) 



~ 2.3 X l{)-^miyE^2Dy^^Jl + z), 



(26) 
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where the last line assumes fiducial plerion parameters. 
We caution that being isotropic in the observer frame, the 
plerion radiation field should actually be highly anisotropic 
in the shock frame moving at F ^ 1. Instead of an ac- 
curate but cumbersome calculation of the Compton up- 
scattered photon distribution including the full angle- 
dependence, we have approximated by using quantities av- 
eraged over angles in the comoving frame; see e.g. Dermer 
1995, Inoue & Takahara 1996, Sikora 1997. 

The broken power-law representations of the EC and 
SSC spectra should be adequate to within an order of 
magnitude, but more accurate calculations properly in- 
tegrating over the broad seed photon frequency distribu- 
tion will result in smoother spectral breaks and somewhat 
larger fluxes (above vf^ and uf^ for fast cooling, or above 



rion at r > i?,,. We fiducially choose n • 



10'''cm"3 



imply- 



,EC 



and 



,sc 



for slow cooling) due to logarithmic terms 



ing that if the; nic;(Uum was globally uniform at this value, 
the blastwave travels from r ~ 3 x lO^^cm to ~ 6 x lO^^cm 
during t ~ 0.2 s to ~ 2 days inside the plerion, illumi- 
nating nearby Fe-rich condensations along the way, and 
turns subrelativistic (F ~ 1) at f ^ 15 days. In reality, a 
sufficiently clumpy baryonic distribution may cause a spo- 
radic type of variability in the afterglow which should be 
detectable by future instruments, and could perhaps be 
relevant to the X-ray flaring behavior observed in the af- 
terglows of GRB970508 (Piro et al. 1999) and GRB970828 
(Yoshida et al. 2001). 

For the case of deceleration occurring mainly within the 
magnetized plerionic material (referred to as the "plerion 
case"), KG02 showed that the usual afterglow formulae 
can be utilized with some modifications, such as substi- 
tuting n with an equivalent baryon density 



1 



4pe + 



{B' + E'f 
2^r 



(27) 



where pe is the electron pressure and B' and E' are 
the magnetic and electric fields in the frame comoving 
with the shocked pulsar wind. For our simplified ple- 
rion model, we estimate rib^eq by taking E' = 0, B' = 

Bp = {2,£,bEp/RIY/'^ and pe = ^eEp/SnRl, which gives 



103cm-3(2^,,o.5/3 + ^B,o.5) t 



p,120- 



contributed by a range of electron Lorentz factors (SEOl). 
Regarding KN effects, a gradual steepening from v <^vkn 
instead of an abrupt cutoff should actually occur, as parts 
of the sec;cl photon spectra will be in the KN regime even 
for electrons with 7 <^ jm- We have not included the 
effects of internal pair attenuation, which could be impor- 
tant for the highest emission energies at relatively early 
times in the afterglow. However, significant differences 
are expected only for TeV energies and above (ZMOl), 
whereby pair degradation by the infrared background dur- 
ing intergalactic propagation should be serious anyway for 
GRBs with ^ ^ 0.1 (e.g. Stecker k de Jager 1998, Totani 
2000). 

3.3. Constraints on the Afterglow Environment 

Having laid out the tools to calculate the broadband af- 
terglow emission, we now discuss some constraints to be 
imposed on the local environment of the blastwave in the 
context of our model. For an afterglow to produce a strong 
Fe line as observed in GRB991216, it is necessary that the 
GRB blastwave decelerates inside the plerion near rela- 
tively dense SNR fragments, so that they are irradiated 
efficiently by afterglow continuum X-rays. 

In case the decelerating medium comprises mostly en- 
trained baryons (referred to as the "baryon case"), we opt 
to directly use the above equations with €3,-2 = 1 and 
ee,-i = 1, since the physics governing these quantities here 
should be similar to that for conventional afterglows de- 
celerating in external media. Although an inhomogeneous 
spatial distribution is actually expected for the baryons, 
we surmise that the formulae are applicable for a certain 
timescale during which the medium can be approximated 
as locally uniform. The local baryon density in the ple- 
rion interior may range anywhere from values as high as 
n ^ lO^^cm"^ for the densest SNR clumps, to much lower 
values n ^ lcm~^ for regions efficiently swept out by the 
magnetized plerionic plasma. If the local density in our 
line of sight n ^ lO^cm"'^, the blastwave turns subrel- 
ativistic in less than ~ Iday and is unable to generate 
the observed afterglows (Lazzati et al. 1999; see however 
Masetti et al. 2001, In 't Zand et al. 2001), whereas if the 
density n ^ 10~^cm~"^, the blastwave does not decelerate 
significantly until it reaches the medium outside the ple- 

'^These conditions may be consistent with the time-delay and ionization parameter constraints discussed e.g. in Piro ct al. (2000) and 
Vietri et al. (2001), since they apply to the distance dj between the continuum source and the irradiated, line-emitting matter, which can 
be smaller than the radius r of the matter from the GRB site. The blastwave (i.e. the continuum source) itself should have moved out to 
r > 10^^ - 10" cm by t ~ Iday, but Fe-rich material could be present in its immediate vicinity at di <g; r. 



Thus, as with 

the baryon case, we fiducially choose n = lO'^cm"^ for 
the plerion case. Spatial uniformity should be a good ap- 
proximation, even for the more sophisticated MHD models 
considered by KG02. One must also make corresponding 
changes for e,. and e^, as well as insert correction factors 
for some quantities, but these depend on the detailed prop- 
erties of the pulsar wind shock (see KG02) for which our 
model lacks information. For want of a better alternative, 
we choose to employ = £,e = 0.5 and £3=^5= 0.5, 
and ignore the correction factors. With these parameters, 
the blastwave should be partially radiative during the fast 
cooling stage; however, an accurate treatment of such cases 
is rather cumbersome (e.g. B5ttcher & Dermer 2000), so 
we proceed with the adiabatic assumption for the sake of 
simplicity. 

A further requirement from the observed Fe lines is that 
tp should be neither much shorter nor much longer than 
^ 120days: ip,i2o ^ 1 would not allow time for enough 
SN-synthesized radioactive Co to decay to Fe (Vietri et al. 
2001), whereas for tj,,i2o ^ 1 the SNR fragments would 
be of too low density to generate the observed Fe line flux, 
unless they can be made extremely condensed (eqs.5,6). 
This justifies our fiducial choice of tp,i2o = 1- ^ 

3.4. Results and Discussion 

Guided by the observed characteristics of afterglows 
such as GRB991216, our fiducial parameters were chosen 
to represent conditions favorable for generating strong Fe 
lines. However, we will not attempt model fits to the avail- 
able data on individual objects, since i) realistically, ad- 
ditional complicating effects of jet geometry, etc. can be 
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important and must be taken into account, and ii) full pa- 
rameter searches for fitting particular observations are out 
of our scope (c.f. Panaitescu & Kumar 2002) The results 
given here for the fiducial parameter sets are intended to 
exemplify some notable features of afterglows wherein EC 
cooling and emission induced by the supranova plerion are 
important. Changing these parameters can lead to either 
stronger or weaker EC effects (see below). 

We first discuss the baryon case. Figs.l and 2 show our 
fiducial results in terms of the time-dependent broadband 
afterglow spectra at different observer times t, and the 
afterglow light curves at selected fixed frequencies, respec- 
tively, assuming the GRB to be at ^ = 1 in the flat lambda 
cosmology (i.e. ~ 2 x lO^^cm). 

The values of X and Y here are fiducially X ^ 1.3 
and Y ~ 2.2, and since they were approximated as be- 
ing constant, the relative luminosities of the synchrotron, 
SSC and EC components stay roughly the same. How- 
ever, each component is prominent in a different energy 
band which changes with t: at the begimiing of the af- 
terglow {t <^ fmin), the SSC component dominates at 
the highest energies TeV) whereas the EC compo- 
nent does so at somewhat lower energies (^ GcV); in 
the latter stages {t ^ Ihr), this is reversed. The lower 
energies are always dominated by the synchrotron com- 
ponent, although the energy of the luminosity peak pro- 
gressively decreases from ~ MeV down to the radio band. 
There are clear differences between the spectral shapes of 
the SSC and EC components: the SSC peaks in vfi, at 
z/'S'C ^ max{i/^^, i^f ^} and has a steep high energy slope 
above this break identical to that of the synchrotron emis- 
sion (energy index — p/2), but the EC component is flat in 



for low frequencies in which tt < < if"-^. 



pfv above ~ max{f, 



EC ,,EC 



}, reflecting the harder 



plerion spectrum. This spectral signature is one way to 
observationally distinguish between SSC and EC emission, 
even though this distinction can be blurred if values of p 

happen to be close to s. 

Here the fast to slow cooling transition occurs at tt ~ 
6.3days, and the afterglow is still fast cooling at f ~ Iday. 
This is longer than in standard discussions, and is partly 
due to the extra cooling process of EC, but primarily 
as a consequence of the relatively high external density 
adopted, n ^ lO'^cm^-*. The light curves at flxed fren- 
quencies can be described by breaks at characteristic times 
tm and tc corresponding to the passage of the minimum 
and cooling frequencies for each spectral component, as 
well as at tt- The synchrotron light curves are as in pre- 
vious studies (Sari et al. 1998): for high frequencies in 
which tc < tjn < tt, U GC t^/^ for t < t^ f„ cx t"^/^ 

for tc < t < tm and OC f(2-3p)/4 [t-l-38] foj. ^ > 

for low frequencies in which tt < tm < tc, fu oc t^/^ for 

t < tt, U OC tl/2 for tt<t< tm, U OC t(3-3f)/4 [^1.13] 
for tm <t <tc and OC i(2-3p)/4 [i-1-38] foj. ^ > 

these expressions and below, the t-dependence inside the 
brackets are for our standard choice of s = 2 and p = 2.5. 
Defining tf*-^, t^ as times when v = ff", i'^ respec- 
tively, the SSC light curves are: for high frequencies in 
which t^'^ < t^ < tt, 



t < t 

iV8 t 



SC 

«^ < i < ti^ 

sc. 

m ' 



(28) 



^1/3 
t 

^(ll-9p)/8 J^- 
^(10-9p)/8 J^- 



t<tt 



J.SC 



1.441 iSC 



1.56 



]t>t, 



<t <t 
SC 



SC 



(29) 



The EC light curves depend on the plerion spectral index 
(-s/2). With tf^, t^ as times when v = v^'^ , v^' re- 
spectively: for high frequencies in which < t^ < U-, 



^1/2 

^(2-3s)/4 [^-1] 
^(2-2p-s)/4 



fECf ^ J.EC 
tg^ <t< tt 



(30) 



'] t > tt; 



for low frequencies in which tt < t^ < tf''^ , 



/.^^oc 



riV2 

^3/2 

f(3-3p)/4 
^(2-2p-s)/4 



^-EC 



t<tt 

tt<t< t^^ 
113] tg^<t<tf^ 
1-2] t > tf^. 



(31) 



Differences can be seen in the high frequency decay in- 
dices at late times between the SSC and EC components: 
/^c OC i(io-9p)/8 both before and after tt for SSC, 

compared to f^^ oc t(2-3s)/4 j^-ij before tt breaking to 
/^<=' a f(2-2p-s)/4 [^-1.2] ^f^gj. foj. EC. These features 
should be important observational diagnostics in the X-ray 
and gamma-ray range, although here again, the discrimi- 
nation is less clear if p approaches s. (Note also that there 
should actually be a slight break in the SSC light curve at 
tt from the time-dependence of Y; SEOl.) The generally 
harder spectrum of the EC emission relative to SSC should 
make it increasingly dominant in the higher energy bands 
at later times. 

The early X-ray and optical emission is predominantly 
synchrotron, but the entry of the minimum frequency of 
the SSC emission (and to a lesser extent, the EC emission) 
gives rise to 'bumps' in the light curves, which appear at 
t ~ 0.5 days in X-rays and t ^ 10 days in the optical 
for our fiducial calculation. In particular, we see that X- 
ray emission after ~ 1 day is SSC plus EC emission. So 
far, there have been no clear indications for such X-ray 
bumps, except for a flattening in the X-ray light curve of 
GRB000926 at t ~ 2 days. Together with multiwavelength 
spectra, this has been interpreted as SSC emission from a 
blastwave in a moderately dense medium of n ~ 30cm~3 
(Harrison et al. 2001), but a different explanation with an 
even denser medium (n ~ 4 x 10**cm~3) is also possible 
(Piro et al. 2001). In the optical band, it is interest- 
ing to recall that a small number of afterglows have been 
seen to exhibit bumps at few tens of days, attributable 
to an underlying supernova light curve superimposed on 
the afterglow (e.g. Bloom et al. 1999, 2002, Galama et 
al. 2000, Lazzati et al. 2001b. The observed bumps seem 
to manifest simultaneous spectral reddening, as opposed 
to the spectral hardening expected for the SSC induced 
bumps. However, with the exception of GRB980326 and 
GRB011121, the combined evidence for light curve bump 
plus reddening is not very strong, so an SSC explanation 
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may not be ruled out. at least for some of the observa- 
tions. In any c:ase, realistie modeling of light curve breaks 
should also include other effects such as jet geometry, so 
firm conclusions in comparison with observations are not 
possible at the moment. 

The detection of the EC component, which is distinc- 
tive to the supranova model, should be best achieved at 
GeV energies and above. We can assess the detectabil- 
ity in GeV gamma-rays by the EGRET and GLAST in- 
struments, following the discussion of ZMOl. The fluence 
threshold of each intrument are described by a constant 
value for photon-limited, short integration times and one 
proportional to t^^"^ for background-limited, long integra- 
tion times. This can be compared with the model fluence 
light curves, intergratcd over the energy range of 400 MeV 
to 200 GeV. Our fiducial baryon case model would have 
been observable by EGRET hours after the GRB from 
z <J 0.5, providing a viable explanation for the extended 
GeV emission seen in GRB940217 (Hurley ct al. 1994). 
Similar observations should also be possible with AGILE 
(e.g. Mereghetti et al. 2001), having a slightly better sen- 
sitivity than EGRET. GLAST should be able to detect 
this emission from redshifts as high as ^; ^ 1.5 up to ^ 1 
day after the burst. 

Our results for the fiducial plerion case are shown in 
Figs. 3 and 4. Since the equivalent baryon density rib^eg = 
lO^cm"^, the only difference here with the fiducial baryon 
case stems from the parameter values of and eb- Many 
of the trends described above are also valid here quali- 
tatively, but important quantitative differences are to be 
seen. The much larger magnetic field (by a factor of 50 
in energy density) induces greatly increased synchrotron 
losses, leading to smaller values of 7c and the correspond- 
ing cooling frequencies, while the larger electron injection 
efficiency (by a factor of 5) causes higher values of 7^ and 
the associated miminum frequencies. This implies very 
large tt, and the electrons remain fast cooling throughout 
the duration of the afterglow. The synchrotron and SSC 
luminosities are larger and occupy higher energy bands; 
the EC emission is relatively much less prominent, appear- 
ing at GeV over the SSC emission only during the early 
time interval from ~ 10 s to ~ 300 s, and again at later 
times <; 10 days. Detection of this component should be 
feasible with future gamma-ray instruments as with the 
baryon case, but here it may be more difficult to entangle 
from the SSC emission. The light curve breaks due to the 
passage of the minimum frequecies occur much later, and 
the SSC bumps are not observable. 

We now discuss the effect of changing our main free 
parameters, particularly those regarding the plerionic en- 
viroment, n and tp. Keeping all others constant for the 
baryon case, decreasing n from n lO^cm"^ leads to 
a greater dominance of the EC component and less syn- 
chrotron emission, as this implies a smaller B and larger X 
for the same cb- Higher values of n may also occur inside 
the plerion; note the recent observations of GRB010222 in- 
dicating n ~ lO^cm"^ for this burst (Masetti et al. 2001, 
In 't Zand et al. 2001), although alternative interpreta- 
tions with lower values of n are possible (e.g. Panaitescu 
& Kumar 2002). More 'conventional' parameters (lower 
n, small X) may be accomodated in the context of our 
model by taking a larger tp to make Up oc Lp/R^ oc tp^ 



smaller, say tp ^ 3yr, and assuming n ^ lcm~^. For 
the plerion case, n is not an independent free parame- 
ter as n = nn,eq oc t~^; here tp ^ 3yr directly implies 
n ~ lcm~^. However, for this larger plerion the volume- 
averaged SN ejecta density should also drop (eq.6), so that 
strong Fe lines are less likely, requiring extreme clumping 
of ejecta material. Note that clear upper limits on Fe line 
emission obtained for some bright afterglows point to a va- 
riety of GRB environments (Yonetoku et al. 2000, 2001), 
which could correspond to a range of tp and/or n in our 
model. We also mention that the well observed afterglow 
of GRB970508 has an Fe line feature detection (Piro et al. 
1999), as well as arguments for consistency with more or 
less 'standard' parameters (Wijers & Galama 1999), but 
including the effects of SSC and EC cooling in the model 
fitting may allow different sets of parameters (e.g. SEOl). 

4. DIRECT DETECTABILITY OF THE PRECURSOR 
PLERION 

Finally, we briefly remark on the possibility of the di- 
rect detection and identification of the precursor plerion 
emission, which would be a straightforward test for the 
supranova model (VS98). A more detailed assessment can 
be found in Guetta & Granot (2002). 

As discussed above, the plerion should be a lumi- 
nous, steady source during a period tp prior to the 
burst, typically with bolo metric luminosity ~ £,eLp ~ 
lO^^erg s~^t~j20 power-law spectrum of energy index 
^ —1 between frequencies ^ lO^^Hz and ^ lO^^Hz. This is 
comparable to high-luminosity quasars and readily observ- 
able; e.g. it should have already been detected out to ^; 
0.7 as bright keV X-ray sources during the ROSAT All- 
Sky Survey (RASS; Voges et al. 1999) The event rate of 
GRBs detectable by BATSE from z ^ 0.7 may be roughly 
Rgrb ~ 30yr-i all-sky (e.g. Porciani & Madau 2001), 
so there should always be at least RcRBtp/fb ^ ^^fb^ 
active plerions all-sky above the RASS sensitivity. Here 
f^^ is inverse of the fractional solid angle subtended by 
the GRB beam, which could be as large as ~ 100 (Frail et 
al. 2001), and we have assumed that the plerions arc not 
strongly beamed. However, the problem lies in ascertain- 
ing these sources as precursors to GRBs, which must be 
shown to be positionally coincident as well as correlated 
in time within tp before the GRB. Even though the op- 
eration periods of BATSE and RASS were close enough 
in time, the large BATSE error boxes {AO ~ 1°) pre- 
clude discrimination from the large number of AGNs with 
similar X-ray spectra. More precise GRB localizations by 
present and next generation instruments (e.g. HETE-2 
with A0 - 10" -10' or Swift with AO - l"-4'), combined 
with concurrently conducted wide field X-ray surveys (e.g. 
XMM-Newton serendipitous survey) should improve the 
prospects on searching for GRB - precursor plerion asso- 
ciations. 

The plerion should also be conspicuous in the optical, 
e.g. its R-band magnitude would be i?c ~ 17 at 2; ~ 1. 
(The precursor supernova may be observable as well, but 
could be masked by the luminous plerion.) Again, identify- 
ing these with the ensuing GRB out of the numerous other 
objects at similar magnitudes requires accurate GRB po- 
sitions and nearly contemporaneous wide field optical sur- 
veys. This may be possible through serendipitous studies 
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of observations made by e.g. the Subaru Suprimc-Cam, 
the Sloan Digital Sky Survey, and the future Supernova 
Acceleration Probe mission. Detection at radio frequen- 
cies may also be possible if there is additional emission 
from the plerion forward shock (§3.1). 

5. SUMMARY AND COMMENTS 

We summarize the salient points of our work. In the 
supranova model of GRBs, the fast-rotating SMNS active 
during the time between the SN and the GRB should drive 
a luminous plerionic nebula into the preburst environment, 
with a number of important consequences for the ensuing 
GRB afterglow. Ray Icigh- Taylor instabilities acting at the 
plerion-SN ejccta interface may induce significant filamen- 
tary clumping of the ejecta material, allowing a favorable 
geometry and enhanced emissivity for the afterglow Fe line 
emission. The plerion radiation field can act as seeds for 
EC cooling and emission in the GRB external shock, lead- 
ing to prominent GeV-TeV gamma-rays during the after- 
glow, which is detectable by GLAST out to typical GRB 
redshifts, and distinguishable from BSC emission by its 
characteristic spectrum and light curve. A direct search 
for the plerion emission prior to the GRB may be con- 
ducted through accurate GRB positions and concurrent 
wide field surveys, e.g. in the optical and X-rays. All of 
these should provide critical tests for the supranova model 
in the near future. 

The ability of the plerion to effectively fragment the 
SNR and penetrate between the ejecta material is an un- 
certain aspect of our scenario. We may therefore consider 
the consequences of the plerion being completely confined 
inside the SNR shell instead. In the case of a spherical 
geometry, such a circumstance would be at odds with con- 
straints derived from the observed Fe lines, since: 1) the 
SNR would likely be optically thick to Thomson scatter- 
ing at the time of the GRB, so that no prompt or after- 
glow emission can be seen (eq.6); 2) the SNR shell, as well 
as the plerionic bubble material, would be dense enough 



to decelerate the GRB blastwave in a time shorter than 
the observed afterglow: and 3) the radiation density of 
the plerion would be so large as to cause excessively fast 
cooling of electrons in the afterglow blastwave. A possi- 
ble alternative is a SNR-confined plerion with a globally 
anisotropic geometry, being elongated in the direction of 
the GRB beam, as discussed by KG02 and Guetta & Gra- 
not (2002). The true efiicacy of RT instabilities at the 
SNR-plerion interface should be clarified through future 
high-resolution numerical simulations including all the rel- 
evant microphysics. Arons (2002) has recently discussed 
in more detail the plausibility of a powerful pulsar wind 
shredding the surrounding SNR in the context of a model 
for cosmic ray acceleration by magnetars. 

Although we have concentrated on the case of the pul- 
sar spinning down exclusively via a magnetospheric wind, 
there is a possibility that significant spindown can also 
occur through gravitational wave emission, such as those 
driven by r-mode instabilities (e.g. Andersson et al. 2000, 
Stergioulas & Font 2001, Lindblom, Tohlinc & Vallisneri 
2001 and references therein). However, many aspects of 
the present theoretical calculations regarding gravitational 
waves are uncertain (Fryer & Woosley 2001). Our work 
can be viewed as exploring the maximum possible effects 
of electromagnetic spindown, which can be modeled with 
some more confidence. In reality, both effects could be im- 
portant to varying degrees, and we leave an investigation 
of such cases to the future. 
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Fig. 1. — The time-dependent afterglow spectra for the fiducial baryon case at t =1 second, 1 minute, 1 hour, 1 day and 15 days after the 
GRB, from top to bottom. The curves for t =1 minute and t =1 day show the decomposition into synchrotron (dotted), SSC (dashed) and 
EC (dot-dashed) components. 
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log t [s] 

Fig. 2. — The afterglow light curves for the fiducial baryon case at i/ = 2.4 x lO^-^Hz (1 GeV, top), ly = 2.4 x lO^'^Hz (1 keV, middle) and 
u = 4.5 X 10^*Hz {Rc band, bottom). The synchrotron (short-dashed), SSC (long-dashed) and EC (dot-dashed) components are also shown 
separately. 
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Fig. 3.— 



As with fig.l, but for the fiducial plerion case. 
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Fig. 4. — As with fig.2, but for the fiducial plerion case. 



